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Abstract 

We are interested in the time asymptotic location of the level sets of solutions to Fisher- 
KPP reaction-diffusion equations with fractional diffusion in periodic media. We show that 
the speed of propagation is exponential in time, with a precise exponent depending on a 
periodic principal eigenvalue, and that it does not depend on the space direction. This is in 
contrast with the Freidlin-Gartner formula for the standard Laplacian. 



Resume 

Propagation dans les equations de type Fisher-KPP avec diffusion fraction- 
naire en milieux periodiques. On s'interesse ici a la localisation asymptotique en temps 
des lignes de niveaux des solutions d'equations de reaction-diffusion de type Fisher-KPP 
avec diffusion fractionnaire en milieu periodique. Nous montrons que la vitesse de propaga- 
tion est exponentielle en temps, avec un exposant dependant d'une valeur propre principale 
periodique, et que cette vitesse ne depend pas de la direction de propagation. Ceci est en 
contraste avec la formule de Freidlin-Gartner pour le laplacien standard. 



Version frangaise abregee 

Considerons l'equation de reaction-diffusion de type Fisher-KPP (Kolmogorov - Petrovskii - 
Piskunov) suivante : 

u t + (-A) Q n = n{x)u -u 2 , x G R d , t > 0, (1) 

avec donnee de Cauchy u(-,0) = uq. Dans ce modele, (— A) a represente le laplacien fractionnaire 
d'ordre a G (0, 1); la fonction \i est supposee periodique en chaque variable X{ d'espace et verifie 
< minu n(x). Les conditions sur la donnee initiale no sont donnees dans le theoreme 
principal. Soit Ai la valeur propre principale periodique de (— A) a — fj,[x)I. Dans le cas Ai ^ 0, 
la solution tend vers quand t — > +oo (voir [3]) ; c'est pourquoi on va supposer Ai < 0. II existe 
alors une unique solution stationnaire strictement positive pour JT]), notee u+. Par unicite, u+ 
est periodique. Si u est solution de ([1]) alors u(x,t) — > u+(x) quand t — > +oo, uniformement sur 
tout compact et on cherche a comprendre a quelle vitesse l'etat stable u+ envahit l'etat instable 
0. 
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Lorsque a, = 1, beaucoup de resultats sont connus. En milieu homogene (i.e. n est constante, 
par exemple = 1 et done ii + = 1), Aronson et Weinberger [I] prouvent que les lignes de 
niveau d'une solution u issue d'une donnee a support compact se propagent asymptotiquement 
en temps a vistesse constante egale a 2, independamment de la direction de propagation. Quand 
u est issue d'une donnee faiblement decroissante Uq, la vitesse asymptotique est exponentielle 
en temps (voir En milieu periodique (i.e. \i est non constante et periodique), si no est a 

support compact, Freidlin et Gartner [10] prouvent, avec des outils probabilistes, qu'il existe une 
vitesse notee w*(e) dans chaque direction de propagation e G S' rf_1 (la definition et l'expression 
de w*(e) sont donnees en (fljl). D'autres demonstrations, utilisant des arguments de systemes 
dynamiques ou d'EDP, sont proposees dans |12) . [8], [2]. 

Lorsque a G (0,1) et dans le cas homogene (u = 1), il est prouve dans [S] , [S| que les 

t 

lignes de niveaux de u se propagent approximativement a la vitesse e d + 2a . La transition entre 
les propagations lineaire (a = 1) et exponentielle (a G (0,1)) est traitee dans [7]. Ici nous 
generalisons ces resultats en milieu periodique. Contrairement au cas du laplacien standard, la 
vitesse de propagation ne depend pas de la direction : 



Theoreme 0.1 Supposons Ai < et considerons u la solution de ([T]) ou la donnee initiale uq 
est continue par morceaux, positive, uq ^ et uq(x) = 0(|xp ( ' c(+2a ^) quand \x\ — > +oo. Alors 
pour tout A G (0,minu) ; il existe une constante c\ > et un temps t\ > (ces deux constantes 
dependent de A et uq ) tels que pour tout t t\ : 

{x e R d | u(x,t) = A} C {x G R d | c A e^* ^ [ar| ^ c^ 1 e^t*}. 

La preuve est basee sur la construction de sous-solutions et sur-solutions explicites. Dans la 
version anglaise de cette note, nous presentons une demonstration complete du cas a < ^ (le cas 
general sera traite dans [1]). L'etape principale est le lemme suivant : 



Lemme 0.2 Supposons Ai < et a < \. Pour des constantes strictement positives a, a, B_, B 
et M , soit 

u ( xt ) = et u(xt) = 

[ } +b_{t)\xf +2a {,) \\ 1 \-l + ^t)\x\ d+2a, 

avec b{t) = (M +B rf+2 ae d+2a t ) — 2a~ et b(t) = (-M |Ai | + B d+2a e d + 2 <* ) — &r . 
II existe une constante M > telle que : 

• si < B < (lAilM -1 )^^ et < a ^ (max^i) _1 (l - BJT^M (Aip 1 ), alors u est sous- 
solution de ([1]) pour t > 0. 

2a 

• pour toiiie constante B > 0, si t = (d + 2a)(2a |Ai |) _1 ln(2 _1 + M |Ai|~ B d+2a ) et a ^ 
(min</>i) -1 (l + 2B d+2a M \\\ \ ), alors u est une sur-solution de (JTJ) pour t > to- 

Un lemme analogue est vrai dans le cas a ^ ^. L'idee sous-jacente a la construction de u et 
u est exposee dans la version anglaise de cette note. 
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1 Introduction: Motivation and main result 



We are interested in the time asymptotic location of the level sets of solutions to the equation 

u t + (-A) a u = n{x)u-u 2 , xeR d ,t>0, (2) 

with initial condition u(-,0) = uq, where a G (0, 1), u is periodic in each Xj-variable and satisfies 
< minu ^ fJ-(x), and ( — A) a is the fractional Laplacian. The nonlinearity fi(x)u — u 2 is often 
referred to as a Fisher-KPP type nonlinearity. 

Let Ai be the principal periodic eigenvalue of the operator ( — A) a — fi(x)I. From [3] we know 
that if Ai ^ 0, every solution to ([2]) starting with a bounded nonnegative initial condition tends 
to as t — > +oo. Thus in this paper we assume Ai < 0. Then, by [3], the solution to (|SJ) tends, as 
t +oo, to the unique bounded positive steady solution to ([2]), denoted by u + . By uniqueness, 

is periodic. The convergence holds on every compact set. Hence, the level sets of u spread 
to infinity for large times, and we wish to understand how fast. To do it, we look for a function 
R e (t) going to +oo as t tends to +oo such that, for every direction e G S^ 1 and every constant 
cG (0,1), 



liminf ( inf u(x,t) I > and limsup I sup u(x,t) = 0. (3) 

t^+co \{x=pe,0^p^R e (ct)} ) t^+oo \{x=pe,p^R e (c' 1 t)} J 

We adopt this slightly unusual definition to cover both linear and exponential propagation. 

The case a = 1 is well studied. In homogeneous media, when the function u is constant (say 
equal to 1), pQ establishes that, if uq is compactly supported, then we may choose R e (t) = It 
regardless of the direction e of propagation. See for the case of slowly decreasing initial 
conditions no (here, essentially, any sufficiently rapidly increasing function becomes an R e (t) for 
some uq). In space periodic media, starting from a compactly supported initial data, Freidlin 
and Gartner [10J have characterised R e (t) by 

c*(e') 

R e (t)=w*(e)t, w*(e)= min — (4) 

where c*(e') is the minimal speed of pulsating travelling fronts in the direction e' . Their proof 
uses probabilistic tools; proofs using dynamical systems or PDE arguments are given in |12| . [8], 

For a G (0, 1) and u = 1 in (J2J) , propagation is exponential in time. Although this fact was 

well noted in physics references, the first mathematically rigorous result is [5], [6], which proved 

t 

the exponential propagation in time; in fact, that R e {t) = e d + 2a in this case. The transition 
between linear propagation (a = 1) and exponential propagation (a G (0, 1)) is examined in [7]- 
See [9] for equations of the type (|2|) with nonsingular integral dispersal. 
In this note we prove the following result : 

Theorem 1.1 Assume that X\ < 0. Let u be the solution to ([2]) with uq piecewise continuous, 
nonnegative, uq ^ 7 and uq(x) = 0(\x\~^ d+2a ^) as \x\ — > oo. Then, for every A G (0, min ji), 
there exist c\ > and a time t\ > (all depending on A and uq) such that, for all t ^ t\, 

{x G R d | u(x, t)=A}c{iGl d c A e^* ^ [x| ^ c^ 1 e^= t }. (5) 
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Theorem 1 1 . 1 1 gives that ([3]) holds with R e (t) = e d + 2a . Thus, spreading does not depend on the 
direction of propagation, and this is in contrast with (|4|) for the standard Laplacian. Moreover, 
the estimate that we obtain is much sharper than that in [5], [6 J for [i = 1. Indeed, to guarantee 
the limits in Q, [6] needed to assume \x\ ^ Ce uit (respectively |x| ^ Ce a2t ) with o\ < d _^ 2 a ^ °2- 
Note that Ai = —1 when fx = 1. When A ^ mmu+, ([5]) can not hold since u(x,t) — ¥ u+{x) as 
t — > +oo. On the other hand one may expect it to hold for A E (0,minu + ). Here we prove it for 
A E (0, min fi); it is easy to see that min fi ^ min u+. 

The proof of Theorem 1 1,1 1 is quite simple: it relies on the construction of explicit subsolutions 
and supersolutions, which are themselves based on a nonlinear transport equation, ([6]), satisfied 
asymptotically by a correctly rescaled version of the solution u. The rest of this note is devoted 
to a full proof in the case a < ^ but the result remain true for all a E (0, 1) and this will be 
explained in [5]. 



2 The proof of Theorem 11.11 

We will from now on assume a < ^. Recall that Ai < denotes the principal periodic 
eigenvalue of the operator (— A) a — [i{x)I and that the corresponding periodic eigenfunction is 
denoted by (f>\. 

Step 1. Let us write u(x,t) = 4>i{x)v{x,t) and define w(y,t) = v(yr(t),t) for r(t) = e d + 2a . Thus, 
for y E R rf and t > 0, w solves 



|Al I z 2qJAjJj 

Wt — ; ; - y ■ w y + e d + 2a 



-A) a w 



Kw 



d + 2a 

where we have used Ai < and we have defined 



<j>x(yr(t)) 



[Ai| w - (pxiyr^w 2 



-<:„ rv I f(r " )) -* ( " ) W»)-.(»))* 



\y - y\ 



-2a\\i\t 



If we formally neglect the term e d + 2a 
we get the transport equation 



(— A) a w — fojyfo)) which should go to as t — > +oo, 



w t 



^-y Wy = \\x\w - <h.(yr(t))v? yeR d ,t>0. 



d + 2a' 

Equation ([6]), completed by an initial datum wq, is solved as: 



(6) 



w(y,t) 



1 



0i(yr(i))|Ai| 1 + e-\W(vi (yr(t))-i-Mvr(t))\\i\ 1 



Taking into account (see for instance [U]) that \x\ d+2a u(x,t) is uniformly bounded from above 



and below (but of course not uniformly in t), it is natural to specialise Wo(y) 



case we have 



74+2^- 111 tmS 



w(y,t) 



Myr(t)) lAip 1 (1 - e-l A il*) + e-l A il* + \y\ 



d+2a ' 
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Since (pi is bounded above and below and t tends to +00, coming back to the function v(x,t) = 
w(xr(i) _1 , t), the idea is to consider the following family of functions modelled by w: 

v(x, t) = —— ° d+2a , u(x, t) = M x )v(x, t), (7) 

|Ai| + b(t) \x\ 

and to adjust a > and b(t) asymptotically proportional to e~' Al '* so that the function u(x,t) 
serves as a subsolution or a supersolution to (J2]). 

Step 2. Thus, let v(x,t) be defined by ([7]). Let us consider the operator K defined by 
Kg(x) = C d , a PV [ MX) ~Jl { ^\ g{x)-~g{x))dx. 

JR d \X — X\ 

To prove the result we need to understand the effect of the operators (— A) a and K on v. From 
simple computations there is a constant D > such that 

|(-A) Q u| ^ Db(t)^v and \Kv\ < Db(t)^v in R d . (8) 

Note that the estimate on K does not hold for a ^ ^ and we must replace the term b(t) d+2a by 

2a — 7 

b(t) d + 2a , where 7 is any positive constant between 2a — 1 and 1. 

We use these estimates to find a and b(t), that we denote by a and b(t), such that u(x,t) = 
4>i(x)v(x,t) is a subsolution to ([2]). Using Ai < 0. note that 

Uu t + (-A) a u- fi(x)u + u 2 = 4>iv t + (f)i(-A) a v - Kv - \Xi\4>iv + 4>jv 2 

< (pi^aT 1 {-ti{t)\x\ d+2a + ^((min^i)- 1 + l)6(t)3^5 (9) 

(lAiP 1 + 6(t) |x| d+2a ) - |Ai| (lAxP 1 + 6(t) |x| d+2Q ) + a^} . 
We define M = D ((min</>i) _1 + l) > and we take a constant > 0. The solution to 

d+2a 

-b'(t) + Mb(t)^ +1 - \Xi\b(t) = 0, 6(0) = (MlAiT 1 + 2a , (10) 

d+2a 

is 6(i) = ^MjAip 1 + B'^e^^ 1 ^ ^ and we have 6(i) < Be~l Al l* < B for all f > 0. 
Thus if we choose I? > and a > such that S < (M" 1 | A.i | ) 2^ and a ^ (max0i) _1 (l — 

2a 1 

B_d+2a m ^ the right hand side of inequality ()10p is less than or equal to 0. 

We make u(x,t) = <fti(x)v(x,t), with a and b(t) denoted by a and b(t), a supersolution in a 
similar fashion. Here we must replace D and M by — D and — M in (fTU|) and (|10p . Given any posi- 

2a 2a | | , d-\-2a 

tive constant B > 0. we take b(t) = (-M |Ai| +B d + 2a e d + 2 « *) — — with initial datum 6(0) = 

(-M jAip 1 +5"^)-^. Defining t = (d + 2a)(2a | A x | ) - 1 ln(2~ 1 +M|Air 1 5^), we 

have£e~ |Al1 * < 6(t) < 2^5 for all t ^ t . Thus, taking a > (min<£i) -1 (l + 2B^ M (Aip 1 ), 
we get a supersolution to ([2|) for t > to- 

/Step 5. Let us prove the main theorem in the case a < \. Due to the assumption on uq for large 
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values of x, we can choose a and B satisfying the conditions obtained in step 2 and such that 
u(x,to) ^ Uq(x). Thus, given any A > and using the maximum principle, we get for t > 

{x G R d \ \x\ > C 2 e^*} C {x G R d \ u(x,t) < A}, 

where C$ +2a = A^B^el^l^omax^. 

The last item to prove is that, for all A G (0,min^), there exist c\ > and a time t\ > 
(both depending on A) such that for t~^t\ 

{x G R d | \x\ < cie^*} C {x G K d | u(a?,t) > A}. (11) 

Since u(-,0) Sj no may not hold, we look for a time ti > max((minii + )~~d~, 2M |Ai| d+2a) such 
that u(x,ti) 2~ 1 minii + for |a?| ^ 1. Moreover, we know that u(x,ti) ^ — — — for 

\x\ ^ 1, where c G (0, 1) is a constant (see the proof of Lemma 2.2 in [6] for the computation). 

d+2a 

Consequently we choose a = — a — and B_ = 2 la d to have u(x, 0) ^ u(x,t\) for all 

2*^ | Ai| max 0i |Ai|tp +1 
x G R d . Taking ti larger if necessary, the requirements for B_ and a in step 2 are satisfied. 
By the maximum principle u(x,t — t\) ^ u(x,t) for t ^ t\. Let us define e = — c d min ^ 1 — and 

max 01 

~d+2a = e -|Ai|ti (Axp 1 ^ -1 , and take x such that |x| < cle^J*. Since b(t) ^ £e~ |Al1 *, for 
t ^ ti we get 

amine/)! amin ^ 

u(x,t) > 1 — ^-rr- = r = £ - 

|Air 1 + : Bel A il*ic'i d+2Q 2 (Aip 1 

Thus, we have found an e > such that (|11|) holds for t > t\ (with c\ and A replaced by c[ 
and e respectively). Note that u(x)u — u 2 ^ (min u)n — ti 2 = ii(min fi — u). Now, we use Lemma 
3.3 in [6] (which also holds for supersolutions), replacing the stable state 1 by min/i. For every 
A G (0, min u), we deduce the existence of a positive constant c\ and a time t\ (both depending 

I Ax | 

on A) such that u(x,t) > A for |x| < cie d + 2a and t ^ t\ . Theorem 11.11 is proved by taking 
c\ = min^i,^" 1 ). 



3 Concluding remarks 

Clearly, the transport equation is a time attractor to the renormalised solution w. In the 
forthcoming paper |3J we will study this question in a more precise way. 

The nonlinearity n(x)u — u 2 is quite special, it yields simple solutions to ©. For more general 
nonlinearities f(x,u), it is more intricate to find sub and supersolutions of the form given by u 
in ([7]), and the proof is less transparent than in the present case. This will be developped in [3J, 
as well as more general diffusions than the fractional diffusion. 
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